Directed transport in a classical lattice with a high-frequency driving 
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We analyze the dynamics of a classical particle in a spatially periodic potential under the influence 
of a periodic in time uniform force. It was shown in [S.Flach, O.Yevtushenko, Y. Zolotaryuk, Phys. 
Rev. Lett. 84, 2358 (2000)] that despite zero average force, directed transport is possible in the 
system. Asymptotic description of this phenomenon for the case of slow driving was developed in 
[X. Leoncini, A. Neishtadt, A. Vasiliev, Phys. Rev. E 79, 026213 (2009)]. Here we consider the 
case of fast driving using the canonical perturbation theory. An asymptotic formula is derived for 
the average drift velocity as a function of the system parameters and the driving law. We show 
that directed transport arises in an effective Hamiltonian that does not possess chaotic dynamics, 
thereby clarifying the relation between chaos and transport in the system. Sufficient conditions for 
transport are derived. 



Transport phenomena in nonlinear systems have at- 
tracted growing interest in the recent decades 
In particular, directed transport in periodic potentials 
U{x + 27r) = U{x) under the influence of an unbiased 
external force E{t) has been subject of research in nu- 
merous papers recently. The studies of the transport are 
motivated, e.g., by its prospects for various technological 
applications There is also much related activity in 
the quantum realm 

Let us concentrate on Hamiltonian systems (no dis- 
sipation). Symmetry analysis allows one to formulate 
necessary conditions of existence of directed motion in 
an ensemble of particles with zero average initial veloc- 
ity 0]. However, what are sufficient conditions to be 
imposed on E{t) and U{x) to guarantee the transport? 
What is the average velocity of transport in an arbitrary 
periodic potential U{x) and force E{t)l Despite a lot of 
research and pubUcations in this field, it seems that it is 
not possible to derive explicit analytical expressions for 
average velocity of a particle as a function of system's 
parameters, apart from the two cases: slow or fast driv- 
ing, where the frequency of perturbation is much smaller 
or faster than the unperturbed frequency of the system, 
respectively. In these cases, it is possible to apply meth- 
ods of classical perturbation theory [ll| . In [l^ , classical 
adiabatic theory was applied to the case of slowly time- 
dependent force. Here we develop classical perturbation 
analysis of the opposite, high-frequency limit, which is 
more interesting from an experimental point of view. We 
derive a general formula for the average velocity of parti- 
cles in a periodic potential under a high-frequency drive 
of an arbitrary form, which provides us with sufficient 
conditions for the transport. 

We consider the system with the Hamiltonian 



p2 

H= — + U{x) - xE{ujt), 



(1) 



i.e. a particle in a spatially periodic potential {U{x + 
27r) = U{x)) influenced by a spatially uniform force 



E{ujt) which is periodic in time and has zero mean: 
E{ujt -f 27r) = E{ujt). We assume the force is chang- 
ing fast: e = i ^ 1. Since the potential is defined up 
to a constant, we fix {U{x))x = 0, where (..)^ means 
averaging over the variable /i. Equations of motion are 
x = P, P = -dU/dx + E{t/e). 

We apply canonical perturbation theory, shifting time- 
dependence to higher order terms in e. The main idea is 
to obtain an effective time- independent Hamiltonian, and 
then to determine how the initial distribution of phase 
points is located in the phase space of this new Hamil- 
tonian. Importantly, the effective Hamiltonian does not 
possess chaotic dynamics. 

We start with several preliminary transformations. In- 
troducing the fast time t — t/e { ~ — e/), the new 
Hamiltonian is (the tilde over the new time is omitted 
from now on) eH = e[-^ -I- U(x) — xElt)]. 

We make a canonical transformation P,x~^p,x using 
a generating function Wi{p, x, t) — x{p + ef{t)). The new 
Hamiltonian is 



cH = eH + e 



dWi 



U{x) 



dt ' 2 

In the following, let us omit the hat over the Hamilto- 
nian, which we denote as an 'intermediate'. We choose 
f{t) ~ {E}, where {..} denotes an integrating operator: 
{F} ^ P^^F{s)ds 
Hamiltonian is 



eH 



illFds),. 



Then, the intermediate 



+ U{x) 



(2) 



We finally make a canonical transformation (x,p) 
{x,p) using the generating function 

W2 — xp+e^S{x,p,t) = xp+e'^Si{x,p,t)+e^S2{x,p,t)+.., 

where all functions Si (to be determined later) are peri- 
odic in time. Variables and the Hamiltonian are trans- 
formed as (^- denotes differentiation over new momen- 
tum p) : 



2 



,dS2 



dSo 



+ 



dp dp 
H{p,x,e) 



P = P 



H{p,x,t,e) + e 



> dSi 
dx 

as 



(3) 



where the new Hamiltonian is eH, and we denote it as 
an 'effective'. We now expand the effective {eH) and the 
intermediate (eH) Hamiltonians in powers of e, and com- 
pare the terms with the same powers of e. This give us 
equations defining the generating function S (see Supple- 
mentary information for details): 



'H{p,x,e) 



+e{Hi{p,x,t) + 



Hoip^x) + e'Hi{p,x) +e^'H2{p,x) + .. 
dS 

H{p,x,t,e) + e— ^ Ha{p,x,t) (4) 
-^) + e\H2{p,x,t) + -^) + .. 



n2 



0, Ui 



Uvl)tU"{x), 



We obtain Hi 
where vi = {/}. 

We see that the effective Hamiltonian coincides with 
the unperturbed one up to the fourth order of the per- 
turbation theory. The expressions for all terms of the 
generating function up to the fourth order can be found in 
Supplementary information. These expressions are use- 
ful for studying various phenomena in periodic potentials 
with high-frequency driving. 

Now, for studying directed transport, we consider an 
ensemble of particles distributed along two unperturbed 
symmetric trajectories with the same energy E = Hq = 
£q (see Fig. [T^). Distribution of the particles is uniform 
in canonical 'angle' variable of the unperturbed Hamilto- 
nian (introducing action-angle variables, constant energy 
implies constant action and the uniform distribution over 
the 'angle' variable seems to be natural). We abruptly 
apply (at a certain initial moment to) the perturbing force 
to this distribution. We shall now derive a formula for 
average velocity of particles in the ensemble, averaging 
over the initial time to as well. Procedure for this is as 
following. After making the canonical transformations 
described above, we obtain the time-independent effec- 
tive Hamiltonian where average velocity of each parti- 
cle is easy to calculate. The transformations depend on 
the initial time to- Because of the asymmetries in the 
function f{t), initial symmetric distribution will not be 
symmetric in the new variables. This will give us the 
transport velocity, nontrivial part of which will survive 
even after averaging over ^o- Quantitatively, we will have 



w+ -t- 



2 

2ir 27r 



V± 



e{E±-Es)sgn[p±{x,£o)] 



uj{E±)uJo{£()) 
(27r)2|p(x,fo)| 



(5) 



dtdx 



Here, v± denote contributions from the unperturbed tra- 
jectories H = So with positive and negative 'old' momen- 
tum, correspondingly, p±{x) are the new (transformed) 



momenta of the particles from these trajectories. E± are 
the new energies of these phase points {E± = W) in the 
effective Hamiltonian, and Eg is the energy of the sepa- 
ratrix of the new Hamiltonian dividing areas of bounded 
and unbounded motion; Q{x) is the Heaviside step func- 
tion (particles coming into the bounded region of phase 
space do not contribute to transport). Then, Co is the fre- 
quency of canonical 'angle' variable in the new Hamilto- 
nian (which gives contribution to the transport velocity), 
while the term ujol\p(x,£o)\ gives us distribution of par- 
ticles in a;— coordinate (recall that we consider uniform 
distribution in 'angle' variable of the unperturbed Hamil- 
tonian). Both frequencies we define using 'old' time, so 
that the drift velocity also is defined using 'old' time. 

For simplicity, we consider initial energies not very 
close to the separatrix, £o ~ £s ^ £• Then the integrand 
of Eq.® simplifies to ± f,g4'gg)V 

Expressions for p,x needed in Eq.® are given by 
Eq.®, with terms up to 4th order being presented in 
Supplementary information. Expanding Eq.® in pow- 
ers of e and using the above-mentioned expressions, we 
get after double integration over x and the initial time, 
in the lowest order in e. 



V = -e^f)tLOo{£o) 



1 d^LOo 1 ^3 



Uo 



2 d£'^ 3 a£:3 



£o 



(6) 



This is one of the main results of this Letter. For the 
particular case of U {x) — ~ cos x we obtain, explicitly. 



-6(f + £o)'^oE^('t) 



48(f + £o)(£o-1)3k5(k) 
+ 6(f + fo)(l + 3£:o^)E2(«;)K(K) (7) 
- 50(15 + 17£:2)E(k)k2(k) + (£o - 1)(3 + h£l)Y:\n) 

where K,E are elliptic integrals of the first and second 
type, correspondingly, and k. — jqij^j- The expression is 
valid for all initial energies not very close to the sepa- 
ratrix (not necessarily large ones). Let us compare this 
result with the earlier results of Q (obtained for the same 
potential V(x) = — cosx, but for large energies). In the 
limit of large energies, we have 



\%£l 



(8) 



For the perturbation used in 5], = i?i cos(wi) + 
i?2 cos(2wt -f a), one has (/'') = — |£'Ji?2sina and we 
have 



45 c-ElE^sma 



16 



p6 



(9) 



which reproduces the result of Q up to the constant fac- 
tor, which is different (note that, in the limit of large 
energies, the distribution which is uniform in canonical 
'angle' becomes uniform in coordinate x). 
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— ^ Integral formula (5) 

• Numerical result 
— *— Asymptotic formula (7) 



FIG. 1: (a) Phase portrait of the unperturbed Hamihonian 
and the initial distribution of particles (two segments of tra- 
jectories drawn in bold thick curves), (b) Drift of the average 
position of ensembles of particles with time in the potential 
U{x) — — cosx under the bichromatic perturbation described 
in the main text. We integrate 100 ensembles of particles, 
each ensemble contains 200 particles being distributed along 
two symmetric unperturbed trajectories with £q = 2.0. Dis- 
tribution is uniform in canonical 'angle' variable of the un- 
perturbed Hamiltonian. To each ensemble, we apply a force 
E{t — ti), where ti is distributed uniformly on (0,27r). Av- 
eraging over all ensembles gives us the average position x{t) 
(for convenience, a normalized position (X) = x/2ne is shown, 
where the normalization is such that the slope of {X){t), with 
time measured in units of fast periods, gives the average veloc- 
ity). The inverse frequency e — 0.03. After a certain transient 
period with beating oscillations of the average positions, the 
asymptotic regime is established. 



Note that in the 'new' (fast) time, the drift velocity 
is of the order of e**. That is why we choose to use the 
fourth-order perturbation theory, even though in the fi- 
nal expressions (|6l7p even vi is not needed. In the lead- 
ing term, the force is entered only via {f^)t, which pro- 
vides a sufficient condition for the nonzero drift velocity 
of the third order of e: {f^)t 7^ (at the energy level 



where 



1 ajwo 

2 a£2 



1 d'uiQ c 
3 dS^ '^0 



transport will be strongly 

supressed, but this condition only influence certain spe- 
cific trajectories). Meanwhile, the arguments based on 
symmetries of the perturbation can only give necessary 
conditions for non- vanishing transport. 

Numerically, we prepare a symmetric distribution of 
initial conditions in the phase space of the unperturbed 
Hamiltonian, and (abruptly) apply the force with an off- 
set ti. Applying the time-dependent force E{t — ti), it is 
important to average over the initial phase ti. Otherwise, 
one can obtain directed transport even in the case of har- 
monic driving [6|] , as a function of the initial phase. In de- 
tail, our procedure is as follows. We prepare N copies of 
a symmetric in coordinate and momentum initial phase- 
space distribution. Specifically, we choose as an initial 
distribution a collection of phase points distributed over 
two symmetric trajectories of the unperturbed Hamilto- 
nian with the energy £ = £q. To each copy, we apply at 
t = a force with an initial offset t,;: Ei{t) ^ E{t — ti), 
where Ei{t) is the force acting on the i— th copy of the 
phase space, and ti is an offset uniformly distributed on 
(0, 27r). We average over all phase points in all N copies 
of the phase space. In other words, we average not only 



FIG. 2: (a)Dependence of the average drift velocity on the 
inverse frequency e for U{x) = — cosx and the perturbation 
described in the main text. Linear fitting give the coefficient 
a = 3.05 for the power-law In |w| = a In e. Deviation from the 
asymptotic value a = 3.0 is partly due to the higher-order 
terms which are not so small in this range of e. (b) Aver- 
age velocity as a function of e^. Squares: integral formula of 
Eq. ([5]). Circles: numerical results. Triangles: asymptotic 
formula of Eq. ([7|. Lines are guides to the eye. Numerical 
results and the integral formula are almost indiscernible at 
this scale. The asymptotic formula deviates from numerics at 
larger values of e, where higher order terms become signifi- 
cant. 



over initial phase-space distribution, but also over initial 
phase of the force. 

For a numerical example, we consider the potential 
U{x) = — cosx and the following perturbation: E{t) = 
-Vi sint - 2V2 sin(2i + 7r/4), where Vi =3, V2 = 0.8. 
This corresponds to 

f{t) ^Vi cost + V2COs{2t + TT / 4), (10) 

and we have: vi{t) = V'l sint + sin(2t + f ), (wi)t = 
^ + ^AP)t = We distribute particles along the 

trajectory of the unperturbed Hamiltonian with £q = 2. 

In Table I, we compare the numerical results for the 
average velocity with predictions of the integral formula 
([5]) and the asymptotic formula ([7]). In FiglTjD, time evo- 
lution of the mean coordinate of ensembles of particles 
is shown for e = 0.03. Slope of the fitted line gives the 
average velocity. In Figl^ti., the dependency of the aver- 
age drift velocity on e is shown. Note that the velocity is 
defined using 'old', original time. In Fig. ^}p, predictions 
of the integral formula ([5|), the asymptotic formula ([7]) 
and numerical calculations are compared. 



TABLE I: Average velocity from rmmerics and theory 



e 


Numerical result 


Integral Eq. ((5| 


Asymptotic Eq. JT]) 


0.01 


-1.004 ■10"'' 


-1.002 ■10-'' 


-9.987 -10-^ 


0.02 


-8.130 10-" 


-8.099 -10-" 


-7.990 10-" 


0.03 


-2.784 -10"^ 


-2.781 -10-^ 


-2.697 -lO"^ 


0.04 


-6.782 -10"^ 


-6.756 -10-^ 


-6.392 -10-^ 


0.05 


-1.372 -lO"" 


-1.364 ■10"'' 


-1.248 lO-* 



To conclude, we derive the formula for average velocity 
of particles in a periodic potential under the influence of 
an unbiased high-frequency force. The average velocity 
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is related to a certain integral momentum of the force, 
thereby the formula is useful for a wide range of perturba- 
tions including (but not limiting to) bichromatic or mul- 
tichromatic harmonic driving, anharmonic driving, etc. 
Moreover, it provides sufficient conditions for directed 
transport, thereby conditions based on symmetry of the 
driving only give necessary conditions of existence of di- 
rected transport. Conceptual aspects of the phenomenon 
are clarified now. Indeed, applying classical perturbation 
theory of the 4-th order, one obtains an effective time- 
independent Hamiltonian. Although the Hamiltonian is 
symmetric in momentum, drift of particles occurs be- 
cause the canonical transformation leading to the new 
Hamiltonian asymmetrically transforms original particle 
distribution into the distribution of particles in the new 
variables. A priori, it was not clear that mechanism of 
transport should be like this. In higher orders of the per- 
turbation theory, one obtains an effective Hamiltonian 
with asymmetry in momentum (13j . It might be pos- 
sible that this asymmetry would be responsible for the 
transport, which is shown not to be the case: generally, 
transport arises within the 4th order of the perturbation 
theory. Importantly, the effective Hamiltonian does not 
possess chaotic dynamics. The full Hamiltonian has only 
a narrow stochastic layer on its phase space in the vicin- 
ity of the separatrix of unperturbed Hamiltonian. This 



layer becomes exponentially small in the high-frequency 
limit and cannot be described in any finite order of the 
perturbation theory. Since we consider energies not very 
close to the separatrix, this chaotic layer is completely 
irrelevant for the dynamics. In other words, chaos is not 
needed for the directed transport. 

Procedure of the averaging used in our work is very 
important in the high-frequency case. Without averaging 
over the initial phase, it is not possible to catch the non- 
trivial part of the transport. That is, considering pertur- 
bation with a fixed initial phase and averaging only over 
the phase space, one gets directed transport even with 
the simple perturbation E{t) = cos(cji — 0o)- So far, 
experiments [l^ only probed transport at fixed initial 
phase of the perturbation Ts'] . We here reveal essential 
features of the nonlinear transport in the high-frequency 
regime that has not been probed in the experiments yet, 
and develop analytical theory for them. The theory can 
be probed in the experiments similar to [3], but with 
a high-frequency perturbation. Our work can therefore 
inspire new experiments in this field. 
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Supplementary information 

We expand the effective Hamiltonian in series in e: 

H{p, X, e) = Ho ip, x) + eHi (p, x) + .. 

ids , 



(11) 



^o(p,.) + ^(e^^+e3_ + ..) + 



+ 



+ 



2 dx^ \ dp dp 
+ enA:p.x)+e—{e^— + e^— + ..) + 



1 d^Hi ( 2dSi gS^a ^2 



2 9a;2 V dp dp 

On the other hand, the effective and the intermediate 
Hamiltonians are related as 



dS 

n{p,x,e) = H{p,x,t,e) + €— 

dS 

= Ho{p,x,t) + e{Hi{p,x,t) + 
+ e\H2{p,x,t) + ^) + .. 

T7 (- ,,^dHof 2dSi sdS2 

= Ho{p,x,t) + —{e — + e' — 



(12) 



dp V dx dx 
ldHl(^dSi ^dS2 \2 



■( 



2 dp^ V dx dx 



■■)■ 



+ ... 



We have 

H = ^ + U{x) + epf{t) + e^i^, Ho = ^ + 
U{x), Hi=pm, H2-i^ 



The term ^^-^^ does not influence the dynamics and 
can be safely omitted. 

Comparing terms of the same order in e, we have (bars 

over p are omitted): 

The zero-order terms in e: 



eO : ■Ho{p,x) = Ho{p,x) = ^ + U{x), (13) 



The first-order terms in e: 



dS 

e': Hi = Hi + -^, Hi=p(/)t=0, (14) 
Si = -pviit), vi{t) = {f}. 



The second-order terms: 



2 ^ ^ dHo dSi „ ^dS2 ^ dHo dSi 
dx dp dt dp dx 

n2 =0, 82 = -V2{t)U'{x), 
V2{t) = {vi} (15) 



The third-order terms: 

?^3 = 0, Ss=jyvs{t)U"{x), vs{t) = {V2} 

The fourth-order terms : 

e4. Ui + 5-^3 ^ 1 927^0 / 9^1 \ 2 

dx dp 2 \ Qp J 

dSj dHodSs dHidSi IdH^ /dSiy 

dt dp dx dp dx 2 dp"^ \ dx J ' 

n, = {-{fv2)t-l{vl}t)U"{x) = l{vl),U"{x). 
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